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ABSTRACT 



The goal of the recently approved space-based LARES mission is to measure 
the general relativistic Lense-Thirring effect in the gravitational field of the spin- 
ning Earth at a repeatedly claimed ~ 1% accuracy by combining its node Q with 
those of the existing LAGEOS and LAGEOS II laser-ranged satellites. In this 
paper we show that, in view of the lower altitude of LARES (h = 1450 km) with 
respect to LAGEOS and LAGEOS II (h ~ 6000 km), the cross-coupling between 
the effect of the atmospheric drag, both neutral and charged, on the inclination 
of LARES and its classical node precession due to the Earth's oblateness may 
induce a 3 — 9% year -1 systematic bias on the total relativistic precession. Since 
its extraction from the data will take about 5 — 10 years, such a perturbing ef- 
fect may degrade the total accuracy of the test, especially in view of the large 
uncertainties in modeling the drag force. 

Subject headings: Experimental tests of gravitational theories; Satellite orbits; 
Spacecraft/atmosphere interactions; Harmonics of the gravity potential field; 
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Introduction 



The LARES (LAser RElativity Satellite) satellite, recently approved_| by the Italian 
Space Agency, should have been launched at the end 2009 with a VEGA rocket in 
a circular orbit in clined by 71.5 deg to the Earth's equator at an altitude otj 1450 km 



( jBarbagallo 



20081 1 . It s goal is a measurement o f the general relativistic gravitomagnetic 



Lense-Thirring effect (ILense and Thirring 



19181 1 due to the Earth's rotation at a repeatedly 



claimed ~ 1% level of accuracy in conjunction with the existing LAGEOS and LAGEOS II 
laser-ranged satellites which fly at much higher altitudes, i.e. h ~ 6000 km. The observable 
is a suitable linear combination of the longitudes of the ascending nodes Q of the three 
satellites, because the gravitomagnetic field of the Earth induces a secular precession on 
such a Keplerian orbital element 
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where G is the Newtonian gravitational constant, L is the Earth's spin angular momentum, 
c is the speed of light in vacuum, a is the satellite's semi-major axis and e is its eccentricity. 
In Tabled] we quote the Lense-Thirring precessions for LAGEOS, LAGEOS II and LARES: 
their magnitudes are of the order of 10 1 — 10 2 milliarcseconds per year (mas yr _1 in the 
following) . 

The much larger classical secular precessions induced on Q by the even zonal harmonic 



^ee on the WEB |http://www.asi.it/SiteEN/MotorSearchFullText.aspx?keyw=LARES 
2 See on the WEB http://www.esa.int/esapub/bulletin/bulletinl35/bull35f_bianchi.pdf 



Actually, the launch date has been postponed to late 2010/early 2011. 
http:/ / www.spacenews.com/civil/100115- asi -expects-budge t-remain-flat-2010.html 



! In its originally proposed configuration ( iCiufolini 



19861 1 the semi-major axis of LARES 



was equal to that of LAGEOS, i.e. a = 12270 km. 
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Table 1: LAGEOS, LAGEOS II and LARES: orbital parameters and node precessions due to 
the terrestrial gr avitomagnetic field and the first two even zonal harmonics for L e = 5.86 x 



10 33 kg m 2 s 1 ( 


M 


cCarthv and Petit 


2004 


), J 2 = 0.00108263538, J 4 = -1.619989 x 10~ 6 . 


(Taplev et al. 


2007 


) 



Satellite a (km) e / (deg) f^T (mas yr 1 ) Qj 2 (mas yr *) Qj 4 (mas yr *) 



LAGEOS 12270 0.0045 109.9 30.7 4.538082658 x 10 8 -2.501490 x 10 5 

LAGEOS II 12163 0.014 52.65 31.5 -8.303252509 x 10 8 9.05051 x 10 4 

LARES 7828 0.0 71.5 118.1 -2.0298207310 x 10 9 2.8925357 x 10 6 



coefficients J^, £ = 2,4,6,... of the multipolar expansion of the terrestrial gr avitational 



19661 ) are a major 



potential accounting for the centrifugal oblateness of our planet (IKaula 
source of systematic uncertainty. They can be written as 

Q oU = J2^iJi, (2) 

1=2 

where the coefficients Clj depend on the Earth's mass M and equatorial radius R, and of 
the orbital geometry of the satellite through a, e and the inclination I of the orbital plane to 
the Earth's equator. Since they have the same temporal signature of the relativistic effect 
of interest, they cannot be subtracted from the signal without affecting the recovery of the 
Lense-Thirring effect itself. Thus, it is of the utmost importance to realistically assess the 
uncertainty in them in order to evaluate their percent impact on the gravitomagnetic shift. 
More specifically, the magnitude of the node secular precessions due to the first two even 
zonals for LAGEOS, LAGEOS II and LARES are listed in Table [H It can be noted that 
the J2— induced rates are of the order of 10 s — 10 9 mas yr" 1 , i.e. seven orders of magnitude 
larger than the Lense-Thirring precessions. The perturbations by J4 are of the order of 
10 5 — 10 6 mas yr -1 , i.e. four orders of magnitude larger than the gravitomagnetic effects. 
Such figures immediately demonstrate the difficulty of determining a smallish relativistic 
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effect with respect to a huge classical one, which needs, thus, to be accounted for with the 
appropriate accuracy (about one part in 10 10 ), or, as combining the three satellite data sets 
aims at (see Section [3]), removed from the signal with the same accuracy. 

Up to now major efforts have been devoted to evaluate the bias due to the lingering 
uncertainty SJ^ in the even zonals according to 



Je 



<^2\n. £ \6J e . 



(3) 



1=2 



A reliable evaluation of such a corrupting effect is made difficult by the fact that the 
relatively low alti tude of LA RES brings into play more even zonals than done by LAGEOS 



and LAGEOS II fllorio 



2009|). 



19871 ) like direct 



Concerning the non- conservative orbital perturbations (IMilani et al. 
solar radiation pressure, Earth's albedo, direct Earth's infrared radiation, atmospheric 
drag, thermal effects like the Yarkovski-Schach and Rubincam ones, they have been so far 
regarded as a minor concern because their dir ect impact on the no de of the LAGEOS-type 



satellites is < 1% of the Lense-Thirring effect (ILucchesi 



2001 



2002). 



In this p aper we want to investigate their indirect effects through the cross-coupling 



(IKaula 



19661 ) 



01 



-Jf 



51 



(4) 



between the zonals-induced node precessions and certain non-gravitational perturbations 
affecting the LARES inclination. We will show that, in particular, the impact of the 
atmospheric drag on /lr may play an important role in the evaluation of the error budget 
of the Lense-Thirring test. Indeed, although the direct secular effect of the atmospheric 
drag on the node vanishes, it is not so for the indirect one due to the non-vanishing secular 
decrease of the inclination which maps onto a node effect. Moreover, we will point out 
that it should not be possible to correct the signal for the measured value of / arc by 
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arc without likely affecting the gravitomagnetic signal of interest itself. Clearly, since the 
non-gravitational perturbations do depend on the particular type of satellite considered 
and since LARES has not yet been launched in orbit, our investigation should not be 
required to be more accurate than it can be in the sense that it must be viewed as a 
reasonable sensitivity analysis pointing out a possible source of potential bias and evaluating 
conservatively the largest possible size of the effect examined. 

The paper is organized as follows. In Section [2] we calculate the secular rate of the 
inclination of a LAGEOS-type satellite induced by a drag force and compute it for LARES. 
In Section [3] we calculate its indirect effect on the Lense-Thirring shift through the secular 
precession due to the even zonal harmonics. We also briefly discuss other non-gravitational 
perturbations which may cause a secular variation of the LARES inclination in Section [4j 
Section [5] is devoted to the conclusions. 



2. The effect of the atmospheric drag on the inclination of LARES 



The Gauss equation for the variation of the inclination I is (IMilani et al.lll987l ) 



dl 
~dt 



r cos u 



iki 



-.A. 



(5) 



where n = a/ GM/a 3 is the un-perturbed Keplerian mean motion, u = g + / is the argument 
of latitude, defined as the sum of the argument of pericentre g, which fixes the position of 
the pericentre with respect to the line of the nodes, and the true anomaly / which reckons 
the instantaneous position of the spacecraft from the pericentre, and A v is the out-of-plane 
component of the perturbing acceleration A. 



The drag force per unit mass is ( King- Held 119871 ) 



A D = --C D ZpVV, 



(6) 
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where 



V = v-V> 



(7) 



is the satellite velocity with respect to the atmosphere; v and V a are the geocentric 
satellite and atmosphere velocities, respectively. The other parameters entering eq. ([6]) are 
the drag coefficient Cp, which depends in a complicated way on the i nteraction between 



the gas of particle s in the surroundings of the satellite and its surface (jAfonso et al. 



1985 



Milani et al. 



19871 ). £ = S/m is the area-to-mass ratiqj of the satellite, and p is the density 



of the atmosphere. 

The velocity of the atmosphere, known as ambient velocity, can be written in terms of 
geocentric inertial quantities as 

V A = cj a x r, (8) 



with 



uja = (1 + W ® = (1 + £Ve k , 



(9) 



where k is the unit vector of the z— axis in an inertial geocentric frame chosen aligned 
with the Earth's angular velocity vector. Note that eq. ([9]) accounts for the fact that the 
atmosphere, in general, does not co-rotate exactly with the Earth; maximum observed 
rom the simplifying assumption of exact co-rotation are of the order of 40% 



deviations 


rom 1 


(Kine-Hele 


1987) 



V 



a = wa {-y i + x 3 , 



(10) 



where i and j are the unit vectors in the reference {xy} plane of the geocentric inertial 
frame which coincides with the Earth's equator; the angle between v, which lies in the 
orbital plane, and Va is the inclination /. 



4 S denotes the spacecraft cross sectional area (perpendicular to the velocity). 
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In order to have the out-of-plane component A v of the drag acceleration evaluated onto 
the un-perturbed Keplerian ellipse, to be inserted into the right-hand-side of eq. (J5]), V 
must be projected onto the n direction of the frame co-moving with the satellite; since 



n = sin I sin Q i — sin / cos Q j + cos / k, 



'111 



then, by choosing f2 = 0, 



Onto the unperturbed orbit 



V a ■ n = —UAXsinl. 



(12) 



x = r cos u cos Q — sin u cos / sin Q, 



(13) 



so that it is possible to obtain 



Va ■ n = -wa 1 " sin / cos u. 



(14) 



Since 



v — v r f + Vt t, 



(15) 



it appears clear that if the Earth's atmosphere did not rotate there would not be any 
out-of-plane component of the drag acceleration which, instead, exists because Va ■ n ^ 
for non-equatorial orbits. Thus, the out-of-plane component of eq. ([6]) is 



A v = — -Cn^pVuAr sin I cos u. 



It turns out that it can be posed (jAbd El-Salam and Senna! 



2004J ) 



(16) 



V = | V — V A I ~ V\/ /cr, /cr ~ 1 



;i7) 



so that 



A v ~ — -Co^pvuj at sin I cos u. 



-9 - 



Concerning the approximations used in eq. (IT7|) and eq. ffT8l) . they are justified since 

fc ^ 1+ (v) 2 - 2 (^) C0S/ ' (19) 

where typically Va — 0.5 km s -1 because of eq. fl5j) and eq. (fTUj) (see also Table ??), and 
v ~ ^jGM/a = 7.1 km s" 1 for orbital heights of about 1400 km. 

By inserting eq. fll8p into eq. fl5]) with the un-perturbed relations 

a(l-e 2 ) 



v = na 



1 + e 2 + 2e cos / 



1 + e cos / ' V 
and integrating over an orbital period Pb by means of 

dt (1 - e 2 f' 2 



(20) 



df, 



(21) 



P b 2tt(1 + ecos/) 2 

one finds that there is a non-vanishing secular rate of the inclination to order zero in the 
eccentricity 

(~jf) = -j C D T,pu A a sin J; (22) 

5 



it agrees with (6.17) byf 



Milani et al 



( 11987I ). In obtaining eq. (1221) we considered the 



atmospheric density p constant over one orbital revolution; since for LARES Pb = 3.7 
h, this is certainly a reasonable assumption. In general, p undergoes many irregular and 



complex variations both in position and time, being large 



by th e heating and cooling of the atmosphere (IKing-Hele 
20M). 



y aff e cted by solar activity and 



1987 



Abd El-Salam and Senna! 



According to Table ??, the inclination of LARES will experience a secular decrease of 



dl\ 



-3 x 10 9 rad yr 1 = —0.6 mas yr 



-i 



(23) 



5 AI in flMilani et al.lll987) is th e shift per revolution; in order to be confronted with eq. 
f l22|) . (6.17) bv Milani et all (119871 ) must be divided by P b = 2n/n. By putting Z 1 and 
v = na one recovers just eq. (l22~l) . 
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Table 2: Relevant physical and orbital parameters of the Earth-LARES system. The quoted 
value for Cd is usually used in lite rature, but it refer s typically to altitudes of some hundreds 



km; at 1450 km it may be larger (IMilani et al. 



19871 ). The value of the area-to-mass ratio E 
has been obtained by using for LARES a diameter of d = 37.6 cm and a mass of m = 400 



kg (http://esamultimedia.esa.int/docs/LEX-EC/CubeSat%20CFP%20issue%201.pdf). The 



value of p is that for the Ajisai satellite (ISengoku et al. 



19961 1 which has a semimajor axis of 



7870 km. Concerning the rotation of the atmosphere, the quoted value has been obtained 
by assuming it is about 20% faster than the Earth itself. 



Parameter 


Value 


Units 


Reference 




3.986004418 x 10 14 


m 3 s~ 2 


fMcCarthv and Petit 2004) 


R® 


6378136.6 


m 


fMcCarthv and Petit 2004) 


■h 


0.00108263538 




CTaDlev et al. 2007) 


«LR 


7828 x 10 3 


m 


fBarbaeallo 2008) 


eLR 







fBarbaerallo 2008) 


-^LR 


71.5 


deg 


fBarbaerallo 2008) 


c D 


2.2 




fAbd El-Salam and Sehnal 2004) 


E 


3 x 10" 4 


m 2 kg -1 


(See caption) 


P 


1 x 10~ 15 


kg m~ 3 


fSensroku et al. 1996) 




8.750538 x 10~ 5 


s^ 1 





where mas stands for milli-arcseconds. Concerning the node, whose Gauss variation 
equation is identical to eq. ([5]) with cosu replaced by sinu/ sin J, it can be shown that there 
are no secular effects induced by the atmospheric drag on it; the first non- vanishing term is 
proportional to e 2 sin 2g. 
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3. The impact of the secular decrease of the inclination on the node 

precession due to the oblateness 

Such a decrease of /lr affects also the secular precession of the spacecraft node due to 
the oblateness of the Earth which is a major corrupting effect for the Lense-Thirring signal. 



Indeed, since 



3 / R\ 2 cos/ / 2 



a bias 



• 3 / R\ 2 sin/ J 2 / dl\ . ,„„. 

occurs. For LARES eq. (j25p yields a shift of 18.8 mas yr -1 over one year; since the Lense- 
Thirring precession of the node of LARES amounts to 118 mas yr" 1 , the cross-coupling of 
the inclination perturbation with the oblateness would yield a systematic error of 16% over 
just one year. 

In fact, the data of LARES will be combined with those of the existing LAGEO S and 



LAGEOS II spacecraft according to the following linear combination of their nodes (jlorio 



2005|) 



^lageos + Ci ^lageos ii + ^lares^ Ci = q.358642219, c 2 = 0.075117522 (26) 

in order to cancel out the impact of the mismodelling £ J 2 and 5J4 of the first two even zonal 
harmonics; the general relativistic prediction of the total Lense-Thirring shift, according 
to the linear combination of eq. fl2"6"|) . is 50.7 mas yr -1 . The com bination of eq. (|2"B"1) is 



based on a strategy put forth for the first time in (lCiufolinilll996l ). It turns out that the 
impact of eq. (I2"2l on eq. ( 12"B"|) is 3% yr -1 . Note that since C\ and c 2 are aimed at removing 
the classical effects by J 2 and J 4 to the needed extent, we released them with nine decimal 
digits, i.e. with the matching accuracy. Indeed, as shown by Table [H the largest classical 
effect is about seven orders of magnitude larger than the Lense-Thirring precessions, and 
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the claimed accuracy of the proposed test is 1%. In obtaining such a result we treated 
the coefficients c\ and c 2 , which depend on the semi-major axes, the eccentricities and 
the inclinations of the three satellites, as constant numbers; let us check if a conservative 

of 



uncertainty of the order 



5Jlr — 1 mas in /lr can affect the numerical values of c\ and 



c 2 at the ninth decimal digit. It turns out that 

|ci(/lr) - ^(/lr + 8h R )\ = 1.7 x 10- 9 , |c 2 (/ LR ) - c 2 (J LR + 5I LR )\ = 3 x 10" 10 . (27) 

Thus, we can conclude that the uncertainty in determining the inclination of LARES is 
not a concern about the accuracy required to compute c\ and c 2 . By repeating the same 
analysis for the semimajor axis olr of LARES, it turns out that an uncertainty 5olr — 1 
cm yields a similar conclusion because 

|ci(a LR ) - C!(a LR + Sa LR )\ = 5 x 1(T 10 , |c 2 (a LR ) - c 2 (a LR + <5a LR )| = 5 x 10 -10 . (28) 

Let us see what could be the impact of the uncertainties in parameters like Co and 
wa on our estimates. For 2 < Co < 2.5 we get a substantially unchanged bias 2.7 — 3.4% 
yr" 1 . By assuming ua = = 7.292115 x 10~ 5 s _1 , i.e. by assuming that the atmosphere 
co-rotates with the Earth, the bias amounts to 2.5% yr _1 . Concerning the approximation of 
eq. (JT7j) used for V, i.e. V = vy/k^ ~ v, it is fully justified in our case. Indeed, (Va/v) cos/ 
appearing in it can be approximated with {ua/u) cos / for a circular orbit; for LARES it 
amounts to 0.03 only, thus yielding v^r = 0.97. It must be noted that the effect of eq. (l2"5j) 
should likely affect the LARES data in full because of the difficulty of realistically modelling 
the drag force, especially Co and p; just to give an idea of the uncertainty in their values 



note that when the solar acti vity is low a typica 



altitude is 2 x 10 kg m ( ISengoku et al 



atmospheric density at about 1500 km 



19961 ) . while for the existing LAGEOS satellite 



6 In fact, it may likely be about one order of magnitude smaller because 5r ~ 1 cm yields 
an uncertainty of approximately 0.3 mas at an height of 1450 km. 
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the drag coefficient is C D ~ 4.9 flAfonso et al.lll985l : Milani et al.lll987f ). The fact that C D 
is larger for LAGEOS than for a lower satellite like LARES is only seemingly contradictory 
(the higher the orbit, the smaller the drag). Indeed, the drag coefficient depends on the 
ratio between the average thermal molecular speed of the atmosphere Vr and the orbital 
speed v of the spacecraft. For relatively low orbits, Vr — 1 km s" 1 is typically smaller 
tharJll v ~ 7.5 km s _1 . On the contrary, at higher altitudes the situation changes because 
v becomes smaller and smaller, while Vt increases rapidly d ue to higher temper atures and 



lower mean molecular weight. In this for LAGEOS, flAfonso et all 

2\ n I /if \ 4 



1985|) 



C 



D 



V T 
v 



2 
15 



Vr 

v 



(29) 



with b ~ 1.4 and ( Vr/v) — 0.8, so that Co — 4. For more details, see the discussion in 



( lMilani et al. 



19871 ). pp.104-107. 



In addition to the neutral part icle drag considere d so far it should also be taken into 



account the charged particle drag flAfonso et al. 



19851 ) due to the fact that a spacecraft 



moving in a gas of electrons and ions tends to acquire an electric charge because of the 
collisions with such particles and also because of the photoelectric effect caused by solar 
radiation. The effect of the charged particle drag can be obtained by re-scaling the one due 
to the neutral particle drag by a multipl icative factor b contain i ng, am ong other things, 



the satellite's potential Vq. According to 



Lucchesi and Paolozzil (l200lf ). it may amount to 



about Vq = —0.3 V for LARES, so that b = 3.1 which implies a 9% yr _1 systematic error 
in the measurement of the Lense-Thirring effect with eq. ( 1261) . It must be pointed out that 
the reduction of the impact of the perturbing accelerations of thermal origin should have 
been re a ched by the LARES team with two concentric spheres. However, as explained by 



Andres! (120071 ) . this solution will increase the floating potential of LARES because of the 



r For LARES (see Tabled] for its orbital parameters) it is v — 7.1 km s 1 . 
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much higher electrical resistivity, so that the evaluations presented here may turn out to be 
optimistic. 

Since the extraction of the relativistic effect would require a multi-year analysis, 
typically At = 5 — 10 yr, the action of the overall atmospheric drag on the LARES 
inclination may be a serious corrupting effect over such timescales. 



Ciufolini et al 



(120091 ) objected that, in fact, the disturbing effect examined would 
not appear in the real data analysis procedure because the inclination along with all the 
other Keplerian orbital elements would bej "measured" arc by arc, so that one should 
only have to correct the signal for the measured value of the inclination; after all, the 
same problems, if not even larger, would occur with the semimajor axes of the LAGEOS 
satellites, which are known to undergo still unexplained secular decrease of 1.1 mm d -1 



(IRubincam 



19821 ) and their consequent mappings onto the node rate s. The problem is that 



while a perturbation Aa pertains the in-plane, radial TZ component (IChristodoulidis et al 



19881 ) of the LAGEOSs' orbits, both the Lense-Thirring node 



precession and the shifts in 



1988|) 



the inclination affect the out-of-plane, normal M component (IChristodoulidis et al. 
of the orbit; thus, even if repeated corrections to the semimajor axis could be applied 
without affecting the gravitomagnetic signal of interest, the same would not hold for the 
inclination. This is particularly true in view of the fact that, for still unexplained reasons 
the Lense-Thirring effect its elf has never been est imated, either as a short-arc or as a 



global parameter. Moreover, 



Ciufolini et al 



( 120091 ) claimed that the recent improvements in 



atmospheric refraction modelling would allow to "measure" the inclination of the LAGEOSs 
satellites at a level of accuracy, on average, of 30 /^as for LAGEOS and 10 /xas for LAGEOS 
II. Firstly, the tracking of a relatively low satellite is always more difficult than for higher 



8 Actually, the Keplerian orbital elements are not directly measurable quantities, contrary 
to, e.g., ranges, range-rates, right ascension, declination. 
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targets, so that caution would be needed in straightforwardly extrapolating results valid for 
LAGEOS to the still non-existing LARES. Second, it is difficult to understand the exact 
sense of such claims because they would imply an accuracy 5r ~ a5I in reconstructing the 
orbits of LAGEOS and LAGEOS II, on average, of 0.2 cm and 0.06 cm, respectively. 



4. Other effects potentially inducing secular variations of the LARES 

inclination 

Among the othe r non- conservative forces acting on the LAGEOS-type satellites, also 



the iRubincaml (jl987l ) effect, due to the anisotropic re-emission of the infrared radiation of 
the Earth along the satellite's spin axis, is important. Such an effect arises from the fact 
that the retroreflectors of LAGEOS have a significant thermal inertia of about 3000 s; since 
such time of thermal response is shorter than the satellite's orbital period = 13526 s and 
larger than the spin period, whose nominal value was about 1 s at the launch epoch, the 
perturbation induces a temperature asymmetry between the hemisphere facing the Earth 
and the one opposite to it. The illuminated hemisphere becomes hotter than the dark one 
only after the Earth has passed its pole, causing a time lag effect accounted for by the 
thermal lag angle 9. A net recoil acceleration directed al ong the satellite spin axis 



occurs. It induces a secular rate of the inclination according to (ILucchesi 



20021) 



In it AR U b is the Rubincam acceleration which depends in a complex way on the physical 
and thermal properties of the satellite and of its array of retro-reflectors, 9 is the thermal 
lag angle, and a z is the component of the satellite's spin along the z axis of a geocentric 
equatorial inertial frame having the x axis along the vernal equinox direction. For LAGEOS 
II the secular inclination rate is of the order of 1.5 mas yr _1 . By assuming for LARES the 
same value of A^ u ^ as for LAGEOS II, i.e. Ar^ ~ —7 x 10~ 12 m s~ 2 , eq. ( 1301 yields an 
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effect of the order of about 0.7 mas yr _1 . In fact, it might be finally smaller because of the 
currently ongoing manufacturing efforts of the LARES team aimed at reducing the impact 
of the non-gravitational pe rturbations of the rmal origin on the new spacecraft with respect 



to the LAGEOS satellites ( iBosco et al. 



20071 ) . Moreover, it will depend on the direction of 



the satellite's spin at the injection in orbit. 



5. Conclusions 

In this paper we have shown that certain subtle non-gravitational perturbations acting 
on the forthcoming LARES satellite may corrupt the claimed goal of performing a ~ 1% 
measurement of the Lense-Thirring effect in the gravitational field of the rotating Earth 
because of the lower altitude of the new spacecraft with respect to the existing LAGEOS 
and LAGEOS II spacecraft. In particular, the interplay between the node precessions due 
to the even zonal harmonics of the geopotential, which are a major source of systematic 
error, and the LARES inclination has been investigated. The atmospheric drag, both in 
its neutral and charged components, will induce a non- negligible secular decrease of the 
inclination of the new spacecraft yielding a correction to the node precession of degree 1 = 1 
which amounts to 3 — 9% yr" 1 of the total gravitomagnetic signal pertaining just the node 
itself. Such a corrupting bias would be very difficult to be modeled. Since the extraction of 
the relativistic signature will require a data analysis of about 5 — 10 yr, the effect examined 
here may yield a degradation of the achievable total accuracy of the test. In principle, also 
the Rubincam effect, of thermal origin, should be taken into account because it can induce 
a non-vanishing secular variation of the inclination. Since both the node and the inclination 
enter the out-of-plane, normal component of the orbit of a satellite, it would not be possible 
to correct for the measured values of the inclinations arc by arc without likely affecting also 
the Lense-Thirring signal itself, especially because it has never been estimated along with 
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the other parameters. 
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